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Abstract 

Variational principles for magnetohydrodynamics were introduced 
by previous authors both in Lagrangian and Eulerian form. In this 
paper we introduce simpler Eulerian variational principles from which 
all the relevant equations of barotropic magnetohydrodynamics can be 
derived. The variational principle is given in terms of six independent 
functions for non-stationary barotropic flows and three independent 
functions for stationary barotropic flows. This is less then the seven 
variables which appear in the standard equations of barotropic mag- 
netohydrodynamics which are the magnetic field B the velocity field v 
and the density p. 

The equations obtained for non-stationary barotropic magnetohy- 
drodynamics resemble the equations of Frenkel, Levich & Stilman [1]. 
The connection between the Hamiltonian formalism introduced in [1] 
and the present Lagrangian formalism (with Eulerian variables) will 
be discussed. 

Finally the relations between barotropic magnetohydrodynamics 
topological constants and the functions of the present formalism will 
be elucidated. 

Keywords: Magnetohydrodynamics, Variational principles 
PACS number(s): 47.65.-ha 
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1 Introduction 



Variational principles for magnetohydrodynamics were introduced by pre- 
vious authors both in Lagrangian and Eulerian form. Sturrock [2j has dis- 
cussed in his book a Lagrangian variational formalism for magnetohydro- 
dynamics. Vladimirov and Moffatt [3] in a series of papers have discussed 
an Eulerian variational principle for incompressible magnetohydrodynam- 
ics. However, their variational principle contained three more functions in 
addition to the seven variables which appear in the standard equations of 
magnetohydrodynamics which are the magnetic field B the velocity field v 
and the density p. Kats [?j has generalized Moffatt's work for compress- 
ible non barotropic flows but without reducing the number of functions and 
the computational load. Moreover, Kats has shown that the variables he 
suggested can be utilized to describe the motion of arbitrary discontinuity 
surfaces [5l |6]. Sakurai [7] has introduced a two function Eulerian varia- 
tional principle for force-free magnetohydrodynamics and used it as a basis 
of a numerical scheme, his method is discussed in a book by Sturrock [2]. A 
method of solving the equations for those two variables was introduced by 
Yang, Sturrock & Antiochos [8]. In this work we will combine the Lagrangian 
of Sturrock ^ with the Lagrangian of Sakurai [7] to obtain an Eulerian 
Lagrangian principle which will depend on only six functions. The varia- 
tional derivative of this Lagrangian will give us all the equations needed 
to describe barotropic magnetohydrodynamics without any additional con- 
straints. The equations obtained resemble the equations of Frenkel, Levich 
&; Stilman [T] (see also [9]). The connection between the Hamiltonian formal- 
ism introduced in [1] and the present Lagrangian formalism (with Eulerian 
variables) will be discussed. Furthermore, we will show that for stationary 
flows three functions will suffice in order to describe a Lagrangian prin- 
ciple for barotropic magnetohydrodynamics. The non-singlevaluedness of 
the functions appearing in the reduced representation of barotropic mag- 
netohydrodynamics will be discussed in particular with connection to the 
topological invariants of magnetic and cross helicities. It will be shown how 
the conservation of cross helicity can be easily generated using the Noether 
theorem and the variables introduced in this paper. 

Due to space limitations this paper is concerned only with barotropic 
magnetohydrodynamics. Variational principles of non barotropic magneto- 
hydrodynamics can be found in the work of Bekenstein &; Oron [lOj in terms 
of 15 functions and V.A. Kats ^ in terms of 20 functions. The authors of 
this paper suspect that this number can be somewhat reduced. Moreover, 
A. V. Kats in a remarkable paper [17j (section IV, E) has shown that there 
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is a large symmetry group (gauge freedom) associated with the choice of 
those functions, this imphes that the number of degrees of freedom can be 
reduced. 

We anticipate apphcations of this study both to hnear and non-hnear 
stabihty analysis of known barotropic magnetohydrodynamic configurations 
im [T2] and for designing efficient numerical schemes for integrating the 
equations of fiuid dynamics and magnetohydrodynamics [l3l HU [151 [E] • 

The plan of this paper is as follows: first we introduce the standard 
notations and equations of barotropic magnetohydrodynamics. Next we 
review the Lagrangian variational principle of barotropic magnetohydrody- 
namics. This is followed by a review of the Eulerian variational principles of 
force-free magnetohydrodynamics. After those introductory sections we will 
present the six function Eulerian variational principles for non-stationary 
magnetohydrodynamics. A derivation of the canonical momenta of the gen- 
eralized coordinates appearing in the Lagrangian allows us to derive the sys- 
tem's Hamiltonian which resembles the Hamiltonian introduced by Frenkel, 
Levich &: Stilman [IJ. This is followed by the derivation of a variational 
principle for stationary magnetohydrodynamics. The discussion related to 
the magnetohydrodynamic topological constants concludes our paper. 



2 The standard formulation of barotropic magne- 
tohydrodynamics 

2.1 Basic equations 

The standard set of equations solved for barotropic magnetohydrodynamics 
are given below: 

^ = Vx{vxB), (1) 
V -B = 0, (2) 

§^ + v-(H = o, (3) 

The following notations are utilized: ^ is the temporal derivative, ^ is 
the temporal material derivative and V has its standard meaning in vector 
calculus. B is the magnetic field vector, v is the velocity field vector and p 
is the fluid density. Finally p{p) is the pressure which we assume depends 
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on the density alone (barotropic case). The justification for those equations 
and the conditions under which they apply can be found in standard books 
on magnetohydro dynamics (see for example [2J). Equation ([1]) describes the 
fact that the magnetic field lines are moving with the fluid elements (" frozen" 
magnetic field lines), equation ^ describes the fact that the magnetic field 
is solenoidal, equation ^ describes the conservation of mass and equation 
([!]) is the Euler equation for a fluid in which both pressure and Lorentz 
magnetic forces apply. The term: 

is the electric current density which is not connected to any mass flow. 
The number of independent variables for which one needs to solve is seven 
{v, B, p) and the number of equations (jll3|4p is also seven. Notice that 
equation ([2]) is a condition on the initial B field and is satisfied automatically 
for any other time due to equation ([1]). Also notice that p{p) is not a variable 
rather it is a given function of p. 

2.2 Lagrangian variational principle of magnetohydrodynam- 
ics 

A Lagrangian variational principle for barotropic magnetohydrodynamics 
has been discussed by a number of authors (see for example [2j) and an 
outline of this approach is given below. Consider the action: 

A = I Cd^xdt, 



^ ^ Pi\v'-eip))-^, (6) 

in which e{p) is the specific internal energy. A variation in any quantity F 
for a fixed position r is denoted as 6F hence: 



6A = J SCd^xdt, 



= ^p{\v^ -w{p)) + pv ■ 5v (7) 

in which w = is the specific enthalpy. 

A change in a position of a fluid element located at a position f at time 
t is given by ^(r, t). A mass conserving variation of p takes the form: 

5p = -V- {pi) (8) 
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and a magnetic flux conserving variation takes the form: 

6B = V x{ix B). (9) 

A change involving both a local variation coupled with a change of element 
position of the quantity F is given by: 



hence 

However, since: 
We obtain: 



AF = + (e • V)F, (10) 
l:i.v = 5v + {i-V)v. (11) 



, dr dAf dE , , 

dt dt dt ^ ' 



^^='%- i^- ^)^=% + i^- V)e - (e • V)^. (13) 

Introducing the result of equations (|8I9|13|) into equation ([7]) and integrating 
by parts we arrive at the result: 

6 A = jd^xpv-iWl 

+ I dt{j dS ■ [-p^lv^ - w{p)) + pv{v ■ a + ^Bx ii X B)] 
- /Af.[-.V.„-^-^-±^MVxB,„, ,14) 

in which a summation convention is assumed. Taking into account the 
continuity equation ([3]) we obtain: 

5 A = jd^xpv-iWl 

+ j dt{j dS ■ [-p^lv^ - w{p)) + pv{v ■ + ^Bx (e X B)] 

^ dv 1 ^ ^ ^ 

d'xi ■ [-pVw -P^- P{v ■ V)v - X (V X B)]}, (15) 

hence we see that if 6 A = for a ^ vanishing at the initial and final times and 
on the surface of the domain but otherwise arbitrary then Euler's equation 
([U is satisfied (taking into account that in the barotropic case = 

Although the variational principle does give us the correct dynamical 
equation for an arbitrary ^, it has the following deficiencies: 
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1. Although ^ is quite arbitrary the variations of 8p and 5B are not. They 
are defined by the conditions given in equation ([8]) and equation ([9]). 
This property is not useful for numerical schemes since ^ must be a 
small quantity. 

2. Only equation ^ is derived from the variational principle the other 
equations that are needed: equation ([T]), equation ^ and equation ^ 
are separate assumptions. Moreover equation ([3]) is needed in order 
to derive Euler's equation (|4]) from the variational principle. All this 
makes the variational principle less useful. 

What is desired is a variational principle from which all equations of 
motion can be derived and for which no assumptions on the variations are 
needed this will be discussed in the following sections. 

The reader should also notice two recent interesting papers by Reinhard 
Prix [20| [2T] which discuss the implications of a time shift r in addition to 
the spatial shift ^(r, t) and also considers the case of multi fluid magneto- 
hydrodynamics . 

3 Sakurai's variational principle of force-free mag- 
netohydrodynamics 

Force-free magnetohydrodynamics is concerned with the case that both the 
pressure and inertial terms in Euler equations (j4]) are physically insignificant. 
Hence the Euler equations can be written in the form: 

'^JlIlJil^JyB^O^ (16) 

47r 

In order to describe force-free fields Sakurai [7j has proposed to represent 
the magnetic field in the following form: 

^ = Vx X V??. (17) 

Hence B is orthogonal both to Vx and V??. A similar representation was 
suggested by Dungey pH', p. 31] but not in the context of variational anal- 
ysis. Frenkel, Levich & Stilman [1] has discussed the validity of the above 
representation and have concluded that for a vector field in the Euclidean 
space 7i? the above presentation does always exist locally but not always 
globally. Also Notice that either x or r/ (or both) can be non single valued 
functions (see [Lj equation 20). 
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Both X ^-iid Tj are Clebsch type comoving scalar fields satisfying the 
equations: 

It can be easily shown that provided that B is in the form given in equation 
()17p . and equation (jlSp is satisfied, then both equation ([1]) and equation 
([2]) are satisfied. Since according to equation ([T6|) both V x B and B are 
parallel it follows that equation (I16p can be written as: 



J . Vx = 0, J . Vr? = 0. (19) 

Sakurai [7] has introduced an action principle from which equation (119]) can 
be derived: 



As = j Csd^xdt, 



B"^ (Vx X Vnf , , 

Taking the variation of equation (j20p we obtain: 

5Cs = ^-(Wxx V?? + Vxx Wr?). (21) 
47r 

Integrating by parts and using the theorem of Gauss one obtains the result: 

/Be B 
dS ■ [{5x^ri - 6r]Vx) x —] + j d^ ■ [([Jx]Vr? - M^x) x — ] 

+ J d^x[6x{^V ■ J) - H^X ■ J)] (22) 

in which / dT, represents an integral along the cut and [5f] represents the 
discontinuity of the variations of non single valued functions. We shall show 
later that x can be defined as a single valued function while r] can be either 
single valued or non single valued. Hence if 6As = for arbitrary variation 
6x, ^f] that vanish on the boundary of the domain (including the cut) one 
recovers the force- free Euler equations (fT9|) . 

Although this approach is better than the one described in equation ([6]) 
in the previous section in the sense that the form of the variations 5x, <5?/ is 
not constrained, it has some limitations as follows: 
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1. Sakurai's approach by design is only meant to deal with force- free 
magnetohydrodynamics; for more general magnetohydrodynamics it 
is not adequate. 



2. Sakurai's action given by equation (|20p contains all the relevant phys- 
ical equations only if the configuration is static {v = 0). If the config- 
uration is not static one needs to supply an additional two equations 
(fTSj) to the variational principle. 



4 Simplified variational principle of non- stationary 
barotropic magnetohydrodynamics 

In the following section we will combine the approaches described in the 
previous sections in order to obtain a variational principle of non-stationary 
barotropic magnetohydrodynamics such that all the relevant barotropic mag- 
netohydrodynamic equations can be derived from using unconstrained vari- 
ations. The approach is based on a method first introduced by Seliger &: 
Whitham [T9|. Consider the action: 



^ = j Cd^xdt, 
C = C1 + C2, 



C. ^ .[| + V.M]-p4 -p^f -£.(VxxV,). (23) 

Obviously Vja, f3 are Lagrange multipliers which were inserted in such a way 
that the variational principle will yield the following equations: 

It is not assumed that z^, a, [5 are single valued. Provided p is not null 
those are just the continuity equation ([3]) and the conditions that Sakurai's 
functions are comoving as in equation (jlSp . Taking the variational derivative 
with respect to B we see that 

i? = S = Vx X Vr/. (25) 
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Hence B is in Sakurai's form and satisfies equation ([2]). By virtue of equa- 
tions (j24p we see tliat B must also satisfy equation ([T|). For the time being 
we have showed that aU the equations of barotropic magnetohydrodynamics 
can be obtained from the above variational principle except Euler's equa- 
tions. We will now show that Euler's equations can be derived from the 
above variational principle as well. Let us take an arbitrary variational 
derivative of the above action with respect to v, this will result in: 

6ffA = j d^xdtp5v-[v-Vv-aVx-P'^'n\+ j dS-5vpv + j d&-5vp[v]. (26) 

The integral / dS ■ Svpv vanishes in many physical scenarios. In the case of 
astrophysical flows this integral will vanish since p = on the flow boundary, 
in the case of a fluid contained in a vessel no flux boundary conditions 
Sv ■ h = are induced (n is a unit vector normal to the boundary). The 
surface integral / on the cut of vanishes in the case that the flow has 
zero cross helicity (see section [7]) since in this case v is single valued and 
[i^] = . In the case that that the flow has non zero cross helicity, ly is 
not single valued (see section [7]), in this case only a Kutta type velocity 
perturbation p]5] in which the velocity perturbation is parallel to the cut 
will cause the cut integral to vanish. 

Provided that the surface integrals do vanish and that 6ffA = for an 
arbitrary velocity perturbation we see that v must have the following form: 

v = {J=Vu + aVx + (3Vr]. (27) 

Let us now take the variational derivative with respect to the density p we 
obtain: 

6pA = j d^xdt6p\^v^ - w - ^ - V ■ Vu] 

+ ids- v5pu + j dt- v5p[v] + j d^xv5p\\l. (28) 



Hence provided that § dS-v5pv vanishes on the boundary of the domain and 
/ • v5p[v] vanishes on the cut of v in the case that v is not single valuecj^ 
and in initial and final times the following equation must be satisfied: 

du 1 _r, , , 

- = -„^-„. (29) 



^Which entails either a Kutta type condition for the velocity or a vanishing density 
perturbation on the cut. 
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Finally we have to calculate the variation with respect to both x and rj this 
will lead us to the following results: 

5^A = J (fxdtSx f^Q^^ + V • {pav) - Vr] ■ J] + j) dS ■ x Vt] - vpa\6x 
+ 1 • X V77 - vpa] [Sx] - J d^xpa6x\ll ' (30) 

6r,A = J d^xdt6r][^^ + V • (pPv) + Vx ■ J] + dS ■ [Vx x ^ - vpP]6ri 
+ 1 . [Vx X ^ - vpP] M - J d^xp(3Srj\ll . (31) 

Provided that the correct temporal and boundary conditions are met with 
respect to the variations 6x and 6r] on the domain boundary and on the cuts 
in the case that some (or all) of the relevant functions are non single valued, 
we obtain the following set of equations: 



da _Vr]- J dp _ Vx • J 
dt p ' dt p ' 

in which the continuity equation ([3j) was taken into account. 



(32) 



4.1 Euler's equations 

We shall now show that a velocity field given by equation (j27p . such that 
the equations for a, /5, x, 1^ satisfy the corresponding equations (I24|29|32|) 
must satisfy Euler's equations. Let us calculate the material derivative of v: 

dv dVv da-L dVx , df5-^ dVrj 
It can be easily shown that: 

dt dt dxk 2 dxk ' 

dVr] ^drj -± drj ^ drj 
V— - V-Ufc- — = -Vvk- 



dt dt dxk dxk 

d^X ^ ^dx _ Y^^^ ^ _Vvfc^ (34) 
dt dt dxk dxk 
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In which is a Cartesian coordinate and a summation convention is as- 
sumed. Equations (j24|29p where used in the above derivation. Inserting the 
result from equations (j34l32p into equation (j33|) yields: 

dt ^^dxk~^ dxk ''^ ^ dxk^ ~^ ^2" 

+ ^((Vt?- J)Vx-(Vx- J)Vr?) 

^1 1 ^ ^ ^ 

= -VvkVk + V(-t;2 - w) + - J X (Vx X V??) 
2 p 

V» 1 - - 

= + -J xB. 35 

P P 

In which we have used both equation ([27]) and equation ([25]) in the above 
derivation. This of course proves that the barotropic Euler equations can be 
derived from the action given in equation ()23p and hence all the equations 
of barotropic magnetohydrodynamics can be derived from the above action 
without restricting the variations in any way except on the relevant bound- 
aries and cuts. The reader should take into account that the topology of 
the magnetohydrodynamic flow is conserved, hence cuts must be introduced 
into the calculation as initial conditions. 



4.2 Simplified action 

The reader of this paper might argue here that the paper is misleading. 
The authors have declared that they are going to present a simplified action 
for barotropic magnetohydrodynamics instead they have added five more 
functions a, (3, x, ^ to the standard set B, v, p. In the following we will 
show that this is not so and the action given in equation (j23p in a form 
suitable for a pedagogic presentation can indeed be simplified. It is easy to 
show that the Lagrangian density appearing in equation (|23|) can be written 
in the form: 

+ ^HB-Bf-{m + ?^ + V-M. (36) 

In which w is a shorthand notation for + aVx + /3Vr/ (see equation (j27p ) 
and B is a shorthand notation for Vx x Vr] (see equation (j25p ). Thus C has 
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four contributions: 



C 

t 



C + C.^ -\- C, -\- ^boundary 

'dv dx 



boundary 



dt 



B) 



+ V - (i^pv). 



(37) 



The only term containing v i^ >C^, it can easily be seen that this term will 
lead, after we nullify the variational derivative with respect to v, to equation 
()27p but will otherwise have no contribution to other variational derivatives. 
Similarly the only term containing B \s Cg and it can easily be seen that this 
term will lead, after we nullify the variational derivative, to equation (I25p 
but will have no contribution to other variational derivatives. Also notice 
that the term Choundary contains only complete partial derivatives and thus 
can not contribute to the equations although it can change the boundary 
conditions. Hence we see that equations (i24l) . equation ([29l) and equations 
([32]) can be derived using the Lagrangian density C[a, j3,Xi'ni^i P\ which 

V replaces v and B replaces B in the relevant equations. Furthermore, 
after integrating the six equations (j24l29l32p we can insert the potentials 
a, /3, Xi Vi ^ iiito equations (f27l) and (f25l) to obtain the physical quantities 

V and B. Hence, the general barotropic magnetohydro dynamic problem 
is reduced from seven equations ()l|3l4p and the additional constraint ([2]) 
to a problem of six first order (in the temporal derivative) unconstrained 
equations. Moreover, the entire set of equations can be derived from the 
Lagrangian density C which is what we were aiming to prove. 



4.3 The inverse problem 

In the previous subsection we have shown that given a set of functions 
a, /3, X, satisfying the set of equations described in the previous subsec- 
tions, one can insert those functions into equation (p7|) and equation ([25|) to 

^ ^boundary also depends on V but being a boundary term is space and time it does not 
contribute to the derived equations 
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Figure 1: A thin tube surrounding a magnetic field line 



obtain the physical velocity v and magnetic field B. In this subsection we 



will address the inverse problem that is, suppose we are given the quantities 
V, B and p how can one calculate the potentials a, /3, Xi Vi The treatment 
in this section will follow closely an analogue treatment for non-magnetic 
fluid dynamics given by Lynden-Bell & Katz |24j . 

Consider a thin tube surrounding a magnetic field line as described in 
figure [H the magnetic flux contained within the tube is: 



in which dl is a length element along the tube. Since the magnetic field lines 
move with the flow by virtue of equation ([T]) both the quantities A$ and 
AM are conserved and since the tube is thin we may define the conserved 
magnetic load: 




(38) 



and the mass contained with the tube is: 




(39) 




(40) 
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in which the above integral is performed along the field line. Obviously the 
parts of the line which go out of the flow to regions in which p = have a null 
contribution to the integral. Notice that A is a single valued function that 
can be measured in principle. Since A is conserved it satisfies the equation: 

By construction surfaces of constant magnetic load move with the flow and 
contain magnetic field lines. Hence the gradient to such surfaces must be 
orthogonal to the field line: 

VA • = 0. (42) 

Now consider an arbitrary comoving point on the magnetic field line and 
denote it by i, and consider an additional comoving point on the magnetic 
field line and denote it by r. The integral: 

f^{r) = ^dl + pii), (43) 

is also a conserved quantity which we may denote following Lynden-Bell &: 
Katz [24| as the magnetic metage. is an arbitrary number which can 
be chosen differently for each magnetic line. By construction: 

Also it is easy to see that by differentiating along the magnetic field line we 
obtain: 

V^jL-B = p. (45) 

Notice that ji will be generally a non single valued function, we will show 
later in this paper that symmetry to translations in p will generate through 
the Noether theorem the conservation of the magnetic cross helicity. 

At this point we have two comoving coordinates of flow, namely A, p 
obviously in a three dimensional flow we also have a third coordinate. How- 
ever, before deflning the third coordinate we will flnd it useful to work not 
directly with A but with a function of A. Now consider the magnetic flux 
within a surface of constant load '&(A) as described in figure [2] (the figure 
was given by Lynden-Bell & Katz [21]). The magnetic fiux is a conserved 
quantity and depends only on the load A of the surrounding surface. Now 
we define the quantity: 
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Figure 2: Surfaces of constant load 



Obviously x satisfies the equations: 



1 = 0, i?^V, 



0, 



(47) 



we will immediately show that this function is identical to Sakurai's func- 
tion defined in equation p7|) . Let us now define an additional comoving 
coordinate r]* since V/i is not orthogonal to the B lines we can choose Vr/* 
to be orthogonal to the B lines and not be in the direction of the Vx lines, 
that is we choose rj* not to depend only on x- Since both Vr]* and Vx are 
orthogonal to B, B must take the form: 



B = AVx X Vt?*. 
However, using equation ([2]) we have: 

V-S = VA- (Vx X V??*) =0. 



(48) 



(49) 



Which implies that j4 is a function of XiV*- Now we can define a new 
comoving function r/ such that: 



ri = I A{x,r]*)dr] 
/o 



di] 



0. 



(50) 
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In terms of this function we recover the Sakurai presentation defined in 
equation (fT7|) : 

^ = Vx X V??. (51) 

Hence we have shown how Xi V can be constructed for a known B, p. Notice 
however, that rj is defined in a non unique way since one can redefine rj for 
example by performing the following transformation: r] ?? + /(x) i^i which 
f{x) is an arbitrary function. The comoving coordinates x^V serve as labels 
of the magnetic field lines. Moreover the magnetic flux can be calculated as: 



$ = j B-dS = J dxdrj. (52) 

In the case that the surface integral is performed inside a load contour we 
obtain: 

HX) = f dxdr^ = xf dr^=[ , ^t;^^^ . (53) 

J A J\ I XK'Imax 'Iminj 

Comparing the above equation with equation (|46|) we derive that r/ can be 
either single valued or not single valued and that its discontinuity across 
its cut in the non single valued case is [r/] = 27r. 

We will now show how the potentials a, (5, v can be derived. Let us 
calculate the vorticity Q of the flow. By taking the curl of equation (|27|) we 
obtain: 

LJ = Vx-u = VaxVx + V^x Vr/. (54) 
The following identities are derived: 

u}-Vx = (V/3 X Vr/) • Vx = -V/? • B, (55) 

c5 . Vr? = (Va X Vx) ■VT] = Va- B. (56) 

Now let us perform integrations along B lines starting from an arbitrary 
point denoted as i to another arbitrary point denoted as r. 

m = - [^^dl + Pii), (57) 
a{r)= — —^dl + a{i). (58) 

Ji B 

The numbers a{i), f3{i) can be chosen in an arbitrary way for each magnetic 
field line. Hence we have derived (in a non-unique way) the values of the 
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a, (3 functions. Finally we can use equation ()27p to derive the function u for 
any point s within the flow: 

u{s) = I {v - aVx - P^rj) ' dr + ^i^), (59) 

J i 

in which i is any arbitrary point within the flow, the result will not depend 
on the trajectory taken in the case that v is single valued. If v is not single 
valued on should introduce a cut which the integration trajectory should 
not cross. 

4.4 Stationary barotropic magnetohydrodynamics 

Stationary flows are a unique phenomena of Eulerian fluid dynamics which 
has no counter part in Lagrangian fluid dynamics. The stationary flow is 
defined by the fact that the physical fields v,B,p do not depend on the 
temporal coordinate. This, however, does not imply that the corresponding 
potentials a, /3, x, are all functions of spatial coordinates alone. More- 
over, it can be shown that choosing the potentials in such a way will lead to 
erroneous results in the sense that the stationary equations of motion can 
not be derived from the Lagrangian density C given in equation (j37p . How- 
ever, this problem can be amended easily as follows. Let us choose a, /5, x, 
to depend on the spatial coordinates alone. Let us choose rj such that: 

T] = fj — t, (60) 

in which 77 is a function of the spatial coordinates. The Lagrangian density 
C given in equation (j37p will take the form: 

t = p{[3- e{p)) - ]^p{Vv + aVx + l3Vf)f - ^(Vx x Vf)f. (61) 

The above functional can be compared with Vladimirov and Moffatt [3] 
equation 6.12 for incompressible flows in which their I is analogue to our [3. 
Notice however, that while (3 is not a conserved quantity / is. 

Varying the Lagrangian L = J Cd?x with respect to v, a, /3, x, p leads 
to the following equations: 

V • (p^) = 0, 
pv-Vx = 0, 
p{ij-Vf] - 1) = 0, 

V ■ ya = , 

P 
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V- V/3 

(3 



P 



1^2 

-V + w. 



(62) 



Calculations similar to the ones done in previous subsections will show that 
those equations lead to the stationary barotropic magnetohydrodynamic 
equations: 

V X {iJ X B) = 0, (63) 



p{v ■ V)v = —Vp{p) + 



(V X B)x B 
47r 



(64) 



5 The Simplified Hamiltonian Formalism 



Let us derive the conjugate momenta of the variables appearing in the La- 
grangian density C defined in equation (j37p . A simple calculation will yield: 



dC 



8 I ^ 
" > dt 



-pa, TTr, 



dC 



^ [ dt 



-pP. (65) 



The rest of the canonical momenta 'Kp^Tia^iTp are null. It thus seems that 
the six functions appearing in the Lagrangian density C can be divided 
to "approximate" conjugate pairs: {v, p),{x-,a),{rj, j3). The Hamiltonian 
density can be now calculated as follows: 



dv 
dt 



dx drj 



P 



1 --.2 

+ —B , 



(66) 



in which v is defined in equation (j27p and B is defined in equation (j25p . 
This Hamiltonian was previously introduced by Prenkel, Levich & Stilman 
[1] using somewhat different variable^. The equations derived from the 
above Hamiltonian density are similar to equations (I24p . equation (j29[) and 
equations (j32p and will not be re-derived here. While Frenkel, Levich & 
Stilman [Ij have postulated the Hamiltonian density appearing in equation 
(|66p . this Hamiltonian is here derived from a Lagrangian. 



'The following notations are used in[T]: A = XiA = 7?,A' = 7r;^,M = 7r^ 
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6 Simplified variational principle of stationary barotropic 
magnetohydrodynamics 

In the previous section we have shown that barotropic magnetohydrodynam- 
ics can be described in terms of six first order differential equations or of an 
action principle from which those equations can be derived. This formalism 
was shown to apply to both stationary and non-stationary magnetohydro- 
dynamics. Although for non-stationary magnetohydrodynamics we do not 
know at present how the number of functions can be further reduced, for sta- 
tionary barotropic magnetohydrodynamics the situation is quite different. 
In fact we will show that for stationary barotropic magnetohydrodynamics 
three functions will suffice. 

Consider equation (|47p . for a stationary flow it takes the form: 



Hence v can take the form: 



v-Vx = 0. (67) 



^ Vy X K 

v = — . (68) 



However, since the velocity field must satisfy the stationary mass conserva- 
tion equation ([3]): 

V • (pv) = 0. (69) 

We see that K must have the form K = \/N, where N is an arbitrary 
function. Thus, v takes the form: 

V = . (70) 

P 

Let us now calculate v x B in which B is given by Sakurai's presentation 
equation (j^5|) : 

Vx X VN 
V X B = ( ) X (Vx X Vt]) 

- >V,xW)^V,. ,71) 

Since the flow is stationary N can be at most a function of the three comov- 
ing coordinates x^P^V defined in subsections 14.31 and 14.41 hence: 

dN-^ dN-^ dN-. 
ViV = — Vx + + — V77. (72) 

ox op or] 
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Inserting equation ([72|) into equation ([7T]) will yield: 



^ I ^ ON ^ 
vxB = -Vx^(Vx X V/i) • Vf?. (73) 
P dp- 



Rearranging terms and using Sakurai's presentation equation (|25p we can 
simplify the above equation and obtain: 

I ^ dN ^ 

vxB = — Vx^(V^--B). (74) 
P on 

However, using equation (j44p this will simplify to the form: 

^ dN 

vxB = -Vx-gj^. (75) 

Now let us consider equation ([T]); for stationary flows this will take the form: 

V X {v X B) = 0. (76) 
Inserting equation ([7^ into equation ([63]) will lead to the equation: 

V(|^) X Vx = 0. (77) 

However, since is at most a function of x, /i, f/ it follows that ^ is some 
function of x- 

dN 

IT, = -^<*)- <™) 

This can be easily integrated to yield: 

N = -^^F{x) + G{x,n)■ (79) 
Inserting this back into equation ()70p will yield: 



(80) 



Vxx (-F(x)V/x + f Vt?) 



P 

Let us now replace the set of variables X; ^ with a new set x\ ^ such that: 

X' = 1 F(x)dx, r?' = (81) 

This will not have any effect on the Sakurai representation given in equation 
(j25]) since: 

= Vx X V?? = Vx X Vf/ = Vx' X Vr/'. (82) 
20 



However, the velocity will have a simpler representation and will take the 
form: 

^ Vx'x V(-^ + G'(x',7?')) 



P 



(83) 



in which G' = p. At this point one should remember that /x was defined in 
equation (j43p up to an arbitrary constant which can vary between magnetic 
field lines. Since the lines are labelled by their x', values it follows that we 
can add an arbitrary function of x'; ^' to /i without effecting its properties. 
Hence we can define a new ^' such that: 

^l' = ^i-G'{^,f)'). (84) 

Notice that ^' can be multi-valued; this will be discussed in somewhat more 
detail in subsection 16. 3i Inserting equation ()84p into equation (j83|) will 
lead to a simplified equation for v: 

,= lK2ll2<L. (85) 

p 

In the following the primes on P, V will be ignored. The above equation is 
analogues to Vladimirov and Moffatt's \3\ equation 7.11 for incompressible 
flows, in which our fj, and x P^^Y the part of their A and ^. It is obvious 
that V satisfies the following set of equations: 

v-Vn = 0, v-Vx = 0, v-Vfi=l, (86) 

to derive the right hand equation we have used both equation and 
equation (j25p . Hence fi,x are both comoving and stationary. As for fj it 
satisfies the same equation as f] defined in equation (|60p as can be seen from 
equation (f62l) . It can be easily seen that if: 



basis = {Vx,yV,yp), (87) 

is a local vector basis at any point in space than their exists a dual basis: 

1 ^ -> ^ -. -> V?7 X Vu B 

dual basis = -CVf] x V/i, V/i x Vx, Vx x Vry) = ( ,v, —). (88) 

P P P 

Such that: 

basisi ■ dual basisj = 5ij, i,j E [1,2,3], (89) 

in which Sij is Kronecker's delta. Hence while the surfaces X) Pj 'H generate 
a local vector basis for space, the physical fields of interest H, B are part of 
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the dual basis. By vector multiplying v and B and using equations ()85|25p 
we obtain: 

vx B = Vx, (90) 

this means that both v and B lie on x surfaces and provide a vector basis 
for this two dimensional surface. The above equation can be compared with 
Vladimirov and Moffatt [3] equation 5.6 for incompressible flows in which 
their J is analogue to our x- 



6.1 The action principle 

In the previous subsection we have shown that if the velocity field v is given 
by equation (j85p and the magnetic field B is given by the Sakurai repre- 
sentation equation (j25p than equations (|l|2|3p are satisfied automatically 
for stationary flows. To complete the set of equations we will show how 
the Euler equations ([4]) can be derived from the action given in equation 
([6]) in which both v and B are given by equation (j85p and equation (j25p 
respectively and the density p is given by equation (l44ll : 



,= V..B = V,.(V,xV,) = |lM^. (91) 

In this case the Lagrangian density of equation ^ will take the form: 

C = 4ift^f-e(,))-!$2^ (92) 

2 p Sir 

and can be seen explicitly to depend on only three functions. Let us make 
arbitrary small variations 5ai = {5x,Sili^t^) of the functions = (x,"?/, /^). 
Let us define the vector: 

C = -^Sa^. (93) 

This will lead to the equation: 

Sa^ = -Vai ■ I (94) 

And by virtue of equation (fTOl) we have: 

Aai = 6a^ + ^ ■ V)ai = 0, (95) 

as one should expect since are comoving with the flow. Making a varia- 
tion of p given in equation ([9T|) with respect to will yield equation ([8]). 
Furthermore, taking the variation of B given by Sakurai's representation 
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([25]) with respect to will yield equation ([9]) . It remains to calculate 5v by 
varying equation (j85p this will yield: 

5v = -—V+-V xipixv). (96) 
P P 

Inserting equations (|8|9|96p into equation ([7]) will yield: 

^ ^ / ^ • V X (f X ^) ^ 1 ^2 
5C = v-V X (p^xv) -^^ - 6p{-v^ + w) 

47r 2 
= v.Wxipixv)-^^^^^^^ + V-{pi){^v^ + w). (97) 

Using the well known vector identity: 

A-V X [C X A) = V ■ {{C X A) X A) + {C X A) X A (98) 
and the theorem of Gauss we can write now equation ([7]) in the form: 

5A = j dt{jdS-[p{ixv)xv-^^^^^^ + {\v^ + w)pi] 

+ J d^xi-[pvxuj + JxB-pV{^v'^ + w)]}. (99) 

The time integration is of course redundant in the above expression. Also 
notice that we have used the current definition equation ([5]) and the vorticity 
definition equation (j54p . Suppose now that SA = for a such that the 
boundary term (including both the boundary of the domain and relevant 
cuts) in the above equation is null but that ^ is otherwise arbitrary, then it 
entails the equation: 

pvxuj + JxB-pV{^v'^ + w)=0. (100) 

Using the well known vector identity: 

Iv(u^) = (u • V)u + iTx (V X w) (101) 

and rearranging terms we recover the stationary Euler equation: 

p{v ■y)v = -Vp+ J X B. (102) 
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6.2 The case of an axi-symmetric magnetic field 

Consider an axi-symmetric magnetic field such that the magnetic field is 
dependent only on the coordinate R which is the distance from the axis 
of symmetry and the coordinate z which is the distance along the axis of 
symmetry from an arbitrary origin on the axis. Thus: 

B = B{R,z). (103) 

Any axi-symmetric magnetic field satisfying equation ([5]) can be represented 
in the form: 

B = VPxV{^) + 27tRB^vA. (104) 

In which (p is the azimuthal angle defined in the conventional way and -B^ 
is the component of B in the (p direction. The function P = P{R, z) is the 
flux through a circle of radius R at height z: 

P{R,z)= [ B -dS = 2tt B,{R',z)R'dR'. (105) 

J{R,z) Jo 

For finite field configurations P will have a maximum Pm = P{Rm, Zm) at 
some Rm,Zm- This circle R = Rm will form a line toroid with the other 
constant P surfaces nearby forming a nested set. There can be several such 
local maxima with local nested toroids in a general configuration but the 
simpler case has just one. 

Let us study the relations between the functions P, B^ and the functions 
X, given in equation ()17p . Assuming that the density p is axi-symmetric 
one can see the magnetic load defined in equation (j40p is also axi-symmetric 
and that the surfaces of constant load are surfaces of revolution around 
the axis of symmetry. From equation ([l6l) we deduce that x = x(-R)^)- 
Expressing equation (I17p in terms of the coordinates i?, z results in: 

B = {-^dzXd^ri)R + {dzXdRV - dRxdzr])4> + {^dRxd^i])z. (106) 

In which dy is a short hand notation for ^ and y is a unit vector perpen- 
dicular to the constant y surface. Comparing equation (jlOOp with equation 
(jl04p we arrive with the set of equations: 

dzP = d,xd427Tr]), BrP = dRxd^{27Tr]) . (107) 

From which we derive the equation: 

dzPdnx - dRPd^x = ^ VP x Vx = 0. (108) 
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Hence: 

P = Fix), (109) 

where F is an arbitrary function. We deduce that P is just another type of 
labelhng of the load surfaces. Thus equation (|107p wih lead to: 

a,(2.,) = ^^^„ = if + (110) 

This should be compared with the result of Young et al. [HI equation 5.1]. 
Substituting the above result in equation ()106p will lead to the equation: 

B<j>(t) = (dzXdRV - dRxdzfi)<i) = Vx x Vt). (Ill) 

This can also be written as: 

B^ = 4>-{Vx^ Vr?) = ((^ X Vx) • V^. (112) 

Hence B^f, is proportional to the gradient of fj along the (j) x Vx direction. 
Since (j) x Vx is known we can integrate along this vector to obtain a non- 
unique solution for fj: 

V = [ a^^dl, (113) 

J |(/>xVx! 

in which dl is a line element along the <j) x Vx line. 

6.3 The case of a magnetic field on a toroid 

Our previous definitions of the surfaces of constant load given in equation 
(|4U|) is ambiguous when the field lines are " surface filling" eg on a toroid and 
give no result when the field lines are "volume filling". At equilibrium B 
and V lie in surfaces (there is an exception when B and v are parallel and fill 
volumes). Our former considerations apply unchanged if these surfaces have 
the topology of cylinders but they need generalization when the surfaces 
have the topology of toroids nested on a line (a similar discussion in which 
non- magnetic fiuids are considered can be found in [22]). We consider a 
surface S spanning that line toroid. Each toroid T will meet S in a loop. 
Consider the magnetic fiux ^(T) through that part of S within the loop 
and the mass enclosed by the toroid m{T). The the mass outside the toroid 
is rh(T) = M — m[T). Now express m as a function m(<^) of the magnetic 
fiux then a definition of magnetic load analogous to that for "cylinders" 
parallel to the axis is: 

dm 
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However, there are now two loads corresponding to the two fluxes associated 
with a given toroid. The other load is obtained by taking a cut across the 
"short" circle section of the torus say of constant 0. The magnetic flux 
through such a cross section may be expressed as a function of the total 
mass m{T) within the toroid and 

A* = (115) 

is a second different load. Of course it is also permissible to reexpress the 
flux ^* as a function of the flux <I> then we find: 

dm _ ^ _ X 
d<s> 

The surfaces of constant A are of course the toroids T which also have A* 
constant. 

A similar problem may arise with the definition of the magnetic metage 
defined in equation (I43p . We may wish to define this quantity using the 
magnetic field B and velocity field v. Since those vectors provide a vector 
basis on the load surface, they can be combined in such a way say: B+^{x)v 
to create a vector which is directed along the large loop of the toroid. (A 
different 7 will leave only twists around the short way.) This combination 
represents an unwinding of the field lines so that they no longer twist around 
the short (long) way. Those loops can be thought as composing the surface 
S. Another surface T,' also composed of such untwisted loops can be so 
chosen that the mass between E and E' and between loads A and A + dA is 
some fixed fraction of ^dX. Such E' form suitable constant metage surfaces 
/i corresponding to partial loads A. Notice that 27r// then describes the angle 
from E turned around the toroid by the short way to reach any chosen point. 
Similar use of the other load A* allows us to define another generalized angle 
/i* measured around the long way. A somewhat less physical approach is 
given below. 

Let us consider a toroid of constant magnetic load. Dungey [181 P- 31] 
has considered the case in which magnetic field lines lie on a torus. He has 
shown that one of the functions (ie rj) involved in the representation (jl7p 
should be non-single valued and therefore a cut should be introduced. 

In order to obtain a simple looking cuts we will replace the previous set 
of functions fi,ri with a new set </>*,?7*, which will be defined as follows: 

, * _ 1-1 + G{x)v 
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Figure 3: A torus of magnetic field lines 

G, fi, / are arbitrary functions of x- Therefore ^ and r/ can be given in terms 
of 0*,ry* as: 

V = 4>* + f{x)v*, 

= -G{x)v + ^{xW, (118) 

7]* can be considered as an angle varying over the small circle of the torus, 
while (p* can be considered as an angle varying over the large circle of the 
torus as in figure ([3]) . On the torus of constant magnetic load the (p* , ij* 
functions have a simple "cut" structure. The above equation can also serve 
as a "definition" of /i. Inserting equations (jllSh into equation (fT7|) and 
equation (j85p will result in the following set of equations: 

B = Vxx V0-+/(x)Vxx Vt,-, 

= aM^^au)^^^. (118) 

P P 

Hence B is partitioned into two vectors circulating along the small and 
large circles of the torus. While v has two vector components one along the 
magnetic field B and another along the small circle. 
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7 Topological constants of motion 

Magnetohydrodynamics is known to have the following two topological con- 
stants of motion; one is the magnetic helicity: 

Km = J B-A(fx, (120) 

which is known to measure the degree of knottiness of lines of the magnetic 
field B |25j . In the above equation A is the vector potential defined implicitly 
by the equation: 

^ = V X i*. (121) 
The other topological constant is the magnetic cross helicity: 

nc= [ B- vd^x, (122) 



characterizing the degree of cross knottiness of the magnetic field and ve- 
locity lines. 

7.1 Representation in terms of the magnet ohydro dynamic 
potentials 

Let us write the topological constants given in equation (jl20p and equa- 
tion ()122p in terms of the magnetohydrodynamic potentials a, /5, x, r], v, p 
introduced in previous sections. 

First let us combine equation p7p with equation (|12ip to obtain the 
equation: 

V X (i'-xVr/) = 0, (123) 
this leads immediately to the result: 

^=XV?? + VC, (124) 

in which Q is some function. Let us now calculate the scalar product B ■ A: 

B ■ A = {Vx^^ri) -VC- (125) 

However, since we have the local vector basis: (VXi Vr/, V^) we can write 
VC as: 

VC = ^Vx+|^VA. + ^Vry. (126) 
ox ofi orj 
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Hence we can write: 

B-A=^iVxxV,)-V, = ^§^. (127) 

Ofj. o[x, y, z) 

Now we can insert equation ()127p into equation (|120p to obtain the ex- 
pression: 

nM = J ^d^idxdr^. (128) 

We can think about the magnetohydrodynamic domain as composed of thin 
closed tubes of magnetic lines each labelled by (x^v)- Performing the inte- 
gration along such a thin tube in the metage direction results in: 

dC 

^d^^ = [C]^,„ (129) 

in which [C]x,»? is the discontinuity of the function Q along its cut. Thus a 
thin tube of magnetic lines in which C, is single valued does not contribute to 
the magnetic helicity integral. Inserting equation (jl29p into equation (|140p 
will result in: 

= j [Qx^vdxdr] = J [C]d<^, (130) 
in which we have used equation (j52p . Hence: 

Kl = ^, (131) 

the discontinuity of C is thus the density of magnetic helicity per unit of 
magnetic flux in a tube. We deduce that the Sakurai representation does 
not entail zero magnetic helicity, rather it is perfectly consistent with non 
zero magnetic helicity as was demonstrated above and in agreement to the 
claims made by Frenkel, Levich & Stilman [T]. Notice however, that the 
topological structure of the magnetohydrodynamic flow constrain the gauge 
freedom which is usually attributed to vector potential A and limits it to 
single valued functions. Moreover, while the choice of A is arbitrary since 
one can add to A an arbitrary gradient of a single valued function which 
may lead to different choices of C the discontinuity value [("] is not arbitrary 
and has a physical meaning given above. 

Let us now introduce the velocity expression given in equation (j27l) and 
calculate the scalar product of B and v, using the same arguments as in the 
previous paragraph will lead to the expression: 

..5.|,V.XV„.V..|||^. (13.) 
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Inserting equation (|132p into equation (|122p will result in: 

nc = J ^d^dxdv. (133) 

We can think about the magnetohydrodynamic domain as composed of thin 
closed tubes of magnetic hnes each labehed by {XiV)- Performing the inte- 
gration along such a thin tube in the metage direction results in: 

y^dfi = (134) 

in which [i^]^,?; is the discontinuity of the function v along its cut. Thus a 
thin tube of magnetic lines in which v is single valued does not contribute 
to the cross helicity integral. Inserting equation (jl34p into equation (jl33p 
will result in: 

= j Hx,vdxdv = I [^]d<^, (135) 
in which we have used equation (j52p . Hence: 

the discontinuity of is thus the density of cross helicity per unit of magnetic 
flux. We deduce that a flow with null cross helicity will have a single valued 
u function alternatively, a non single valued v will entail a non zero cross 
helicity. Furthermore, from equation (j29p it is obvious that: 

M - 0. (i3r, 

We conclude that not only is the magnetic cross helicity conserved as an 
integral quantity of the entire magnetohydrodynamic domain but also the 
(local) density of cross helicity per unit of magnetic flux is a conserved 
quantity as well. 

In the following sub section we give a simple example which will demon- 
strate some of the general assertions of this paragraph. 

7.2 A Helical Stratified Magnetic Field 

Consider a magnetohydrodynamic flow of uniform density p. Furthermore 
assume that the flow contains a helical stratified magnetic field: 

-^i2B^{l-^)4> + B,oz R<a , . 

I R> a ^ ^ 
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in which R,(l),z are the standard cyUndrical coordinates, R,<j),z are the 
corresponding unit vectors and Bzo,B± are constants. The magnetic field 
is contained in a cyhnder of Radius a and is independent of z. A possible 
choice of the vector potential A is: 

X=lB.oxy + BMl-^r^, R<a 

10 R> a ^ ' 

in which y is a unit vector in the y direction. Let us calculate the magnetic 
helicity of the field using equation (jl20p . In order to obtain a finite magnetic 
helicity we assume that the field is contained between the planes z = and 
z = 1 furthermore we assume that the planes z = and z = 1 can be 
identified such that the magnetic field lines are closed. Thus the domain 
becomes a topological torus. Inserting equation ()138|) and equation ()139p 
into equation (|120|) will result in: 

Hm = J B- Ad^x = |a35,o5± (140) 

First let us calculate to load using equation (|10]) (we assume that R < a in 
the following calculations) we obtain: 

hence the load surfaces are cylinders. The x function can now be calculated 
according to equation ([16|) to yield the value: 

X = \b,oR\ (142) 
Solving equation (jSip for r] we obtain the following non unique solution: 

r/=^(l--)-+<^. (143) 
BzQ a a 

Substituting equation (jl42p , equation (jl43p and equation (|139p into equation 
(jl24p we can solve for Q and obtain: 

C = B^z{a -R) + ^B,oxy. (144) 

Since we have identified the z = and z = 1 planes the z coordinate is not 
single valued and therefore (" is a non single valued function which has a 
discontinuity value: 

[(] = B^{a-R). (145) 
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Thus we can calculate the magnetic helicity using equation ()130p and obtain: 

= y [C]d^ = ^a^B,oB^, (146) 
which coincides with the result of equation p4U|) . 

7.3 Cross Helicity Conservation via the Noether Theorem 

The conservation of helicity / v • Cod'^x in ideal (non-magnetic) barotropic 
fluid when certain conditions are satisfied in particular, when a; • n = on 
the (Lagrangian) surface bounding the volume of integration was discovered 
by Moffat [25] . Moreau [26] has discussed the conservation of helicity from 
the group theoretical point of view. In his paper he used an enlarged Arnold 
symmetry group [27] of fluid element labelling to generate the conservation 
of helicity. Yahalom [29\ [30j has shown that the symmetry group generating 
conservation of helicity becomes a very simple one parameter translation 
group in the space of labels (alpha space) when represented by Lynden-Bell 
and Katz [24] labelling. We will now show that in the case of magnetohy- 
drodynamics the same one parameter translation group will generate the 
magnetic cross helicity via the Noether Theorem. 

Let us denote the initial position of a fluid element by (xq), by mass 
conservation: 

pix')d'x = p{x^o)d'xo = p(xo^)§444l^'«- (147) 

Since the initial position of a fluid element can not depend on time it must 
depend on the label only, and therefore by an appropriate choice of the a' s 
we obtain: 

p{x'^)d^x = d^a, p = ^) / ' „( . 148 

Where we assume that the above expressions of course exist. Let us look at 
the action A defined in equation ([6]). 

Prom the discussion following equation (llSp we know that if the ^ vari- 
ations disappear at times tQ,ti than A is extremal only if Euler's equations 
are satisfied and the boundary term disappears. If on the other hand we 
make a symmetry displacement i.e. a displacement that makes 5A van- 
ish and assume that Euler's equations are satisfied and the boundary term 
disappears, we obtain that: 

/ V ■ id^a = const. (149) 
JV 
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this is Noether's theorem in it's fluid mechanical form. 

The a' s so chosen as to satisfy the equation (jl48p are not unique in fact 
one can always choose another set of variables say a' s such that: 

It is quite clear that if the domain of integration is not modified any new set 
of q's satisfying equation (jl50p can be chosen without effecting the value of 
the Lagrangian L. This is nothing but Arnolds |27jj alpha space symmetry 
group under which L is invariant (see also Katz &; Lynden-Bell [28j). For 
some flows the domain of integration can be modified without effecting L, in 
that case we have additional elements in our symmetry group. If we make 
only small changes 6a than we can define the group as follows: 



U, da ■ n 



= 0, (151) 

surface 



dak 

where n is a unit vector orthogonal to the surface of the alpha space volume 
which we integrate over. The restriction 5d ■ n = is only needed 

surface 

when the infinitesimal transformation changes the domain of integration in 
such a way as to modify L. In this paper we are interested in the subgroup 
of translation, i.e.: 

5ak = flfc, flfc = const. (152) 
This subgroup of course does not satisfy da ■ n =0 unless at least 

surface 

few of the a's are cyclic or L is not effected by the modification of domain. 

In section l4?3l we have defined the three following parameters: the mag- 
netic load A, the magnetic metage fi and ij. Notice that since the magnetic 
lines are closed fi is an angular variable and we can translate it with out 
changing L. Choosing a'^ = X) and inserting those variables into equa- 
tion (jl48p we re-derive equation ([9T|) . 

The appropriate ^ symmetry displacement associated with the infinites- 
imal change in a^ is given by equation (j93p . For a metage displacement ^ 
takes the form: 

i=-^Sti = -5ti-. (153) 

Inserting this expression into the boundary term in equation (jl5p will result 
in: 

SAb = j dtjdS- [B{]^v'^ - w{p)) - v{v ■ B)] = 0, (154) 
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which is indeed Moffat's condition for magnetic cross hehcity conservation 
|25j as expected. Inserting equation (jl53p into equation (jl49p we obtain the 
conservation law: 



Thus we conclude that the alpha translation group in the direction of // gen- 
erates conservation of helicity. [ One could of course introduce the symmetry 
displacement ^ = e-^, however, in this case one should show that the above 
displacement is a symmetry group displacement which is not obvious if we 
do not take into account Arnold's group and Lynden-Bell & Katz labelling. 
Moreover in coordinate space the symmetry group appears arbitrary com- 
plex depending on the flow considered as opposed to its apparent simplicity 
in alpha space.] 

8 Conclusion 

In this paper we have reviewed variational principles for barotropic magneto- 
hydrodynamics given by previous authors both in Lagrangian and Eulerian 
form. Furthermore, we introduced our own Eulerian variational principles 
from which all the relevant equations of barotropic magnetohydrodynamics 
can be derived and which are in some sense simpler than those considered 
earlier. The variational principle was given in terms of six independent 
functions for non-stationary flows and three independent functions for sta- 
tionary flows. This is less then the seven variables which appear in the 
standard equations of magnetohydrodynamics which are the magnetic field 
B the velocity field v and the density p. 

The equations in the non-stationary case appear have some resemblance 
to the equations deduced in a previous paper by Frenkel, Levich and Stil- 
man [1]. However, in this previous work the equations were deduced from 
a postulated Hamiltonian. In the current work we show how this Hamil- 
tonian can be obtained from our simplified Lagrangian using the canonical 
Hamiltonian formalism. 

The appearance of a non zero magnetic helicity and cross helicity, is 
connected with the fact that some of the functions which we defined are 
non-single valued. This was elaborated to some extent in the final section of 
this paper and was connected to the properties of the functions C, i^- We have 
also shown that the density of cross helicity per unit of magnetic flux is also 
a conserved quantity and is equal to the discontinuity of v. Furthermore, we 




(155) 
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have shown that the conservation of cross hehcity can be deduced using the 
Noether theorem from the symmetry group of magnetic metage translations. 

The problem of stability analysis and the description of numerical schemes 
using the described variational principles exceed the scope of this paper. We 
suspect that for achieving this we will need to add additional constants of 
motion constraints to the action as was done by [311 [32] see also [33] , hope- 
fully this will be discussed in a future paper. 
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